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We determine the pseudo-critical couplings at imaginary chemical potentials by high-statistics 
Monte Carlo simulations of QCD with four degenerate quarks at non-zero temperature and baryon 
density by the method of analytic continuationan. We reveal deviations from the simple quadratic 
dependence on the chemical potential visible in earlier works on the same subject. Finally, we 
discuss the implications of our findings for the shape of the pseudo-critical Une at real chemical 
potential, comparing different possible extrapolations. 
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1. Introduction 

Understanding the phase diagram of QCD on the temperature-chemical potential (r,/x) has 
many important implications in cosmology, in astrophysics and in the phenomenology of heavy ion 
collisions. Unfortunately, the study of QCD at nonzero baryonic density by numerical simulations 
on a space-time lattice is plagued by the well-known sign problem: the fermion determinant is 
complex and the Monte Carlo sampling becomes unfeasible [1]. 

One of the possibilities to circumvent the sign problem is to perform Monte Carlo numerical 
simulations for imaginary values of the chemical potential, where the fermion determinant is real 
and the sign problem is absent, and to infer the behavior at real chemical potential by analytic con- 
tinuation. The idea of formulating a theory at imaginary chemical potential ix was first suggested 
in Ref. [2]. Soon after there were first appUcations to QCD [3-7]. 

The method of analytic continuation has several advantages, indeed coupling j3 and chemical 
potential n can be varied independently and there is no Umitation for increasing lattice sizes. But 
the extent of the attainable domain with real /i is limited by the periodicity and non-analyticities [8] 
and by the accuracy of the interpolation of data for imaginary jj.. 

Analytic continuation can be exploited for evaluating physical observables at real chemical 
potential. Indeed a careful numerical analysis in SU(2) has shown that a considerable improvement 
can be achieved if ratio of polynomials are used as interpolating function [9]. However the main 
goal of the application of the method of analytic continuation is locating the critical fine on the 
(r,/i)-plane for real jj.. 

In previous studies [10, 11] we investigated two-colors QCD at finite baryon chemical poten- 
tial and QCD at finite isospin chemical potential. These theories are free from the sign problem 
and simulations can be performed both at real and imaginary chemical potential. In this way it 
is possible to check the reliability of the analytic continuation. The lesson we learned from the 
aforementioned studies was that the prediction for the pseudocritical couplings at real chemical 
potentials may be wrong if data at imaginary n are fitted according to a linear dependence. 

We present here results obtained in the determination of the pseudocritical line in SU(3) with 
Nf = 4a.t finite baryon density [12]. 

2. Numerical results 

We have considered consider QCD with Nf = 4 degenerate standard staggered fermions of 
mass am = 0.05. Lattice simulations have been performed on a 12^ x 4 lattice using the exact ^ 
algorithm [13], properly modified for the inclusion of a finite chemical potential. Typical statistics 
collected is 10k trajectories of 1 molecular dynamics unit, growing up to 100k trajectories for the 
values of the couplings around the peak of the susceptibility of a given observable. 

In order to obtain the critical line on the (r,;U)-plane for real ^ we first need to locate the 
(pseudo-)critical couplings for several fixed values of the imaginary chemical potential. This has 
been done by looking for peaks in the susceptibilities of a given observable. A suitable interpola- 
tion of the (pseudo-)critical couplings at imaginary chemical potential is then extrapolated to real 
chemical potential. 



2 



Critical line of QCD 



L. Cosmai 



In SU(3) with Nf = 4 the critical line is a line of first order transitions in the first Roberge- 
Weiss sector —(7r/3)^ < (/i/T)^ <0. Therefore tunneling between the different phases is expected 
every few thousands trajectories (Fig. |l|). 
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Figure 1: (Left) Distribution of the real part of the Polyakov loop in SU(3) with Nf =4 on a 12^ x 4 lattice 
with am=0.05 at ajj, = 0.170/ and for two j3 values around the transition. (Right) Monte Carlo history of the 
real part of the Polyakov loop in SU(3) with Nf — 4 ona 12^ x 4 lattice with a;«=0.05 at an = 0.170i and 
J3=5.066. 




Figure 2: Susceptibility of the (real part of the) Polyakov loop vs j3 in SU(3) with A'/ = 4 on a 12^ x 4 lattice 
with am=0.05 and ajJ. — 0.150/. The solid lines represent the Lorentzian interpolation. 



The critical coupling jS(/i^) at imaginary chemical potential /i is determined as the value for 
which the susceptibility of (the real part of) the Polyakov loop exhibits a peak (Fig. ^. The peak 
value at given imaginary chemical potential has been determined through a Lorentzian fit. We 
also checked the consistency of our determinations by means of Ferrenberg-Swendsen method or 
by estimating the point where the peaks in the distribution of the real part of the Polyakov loop 
has equal height. Moreover we verified the determination of the pseudo-critical couplings do not 
change if other observables are used to perform the aforementioned analyses. 

In Fig. ^the values of the critical couplings j8(/x^) versus ajJ.^ are displayed. Data do not line 
up along a straight line and ^{ix^) cannot be parametrized by a polynomial of order jx^. Therefore 
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Figure 3: Critical couplings obtained in SU(3) with A'/ = 4 on a 12^ x 4 lattice with flm=0.05. The dashed 
vertical line indicates the boundary of the first RW sector, fllm(/i) = n/l2. 



we are lead to use the following fit function: 



(2.1) 



We found that a ;^^/dof ~ 1 is obtained both with sextic polynomial in pL (Fig. 4 (left)) and ratio of 
polynomials 2nd to 4th order or 4th to 2nd order (Fig. 4 (center)). 
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Figure 4: Fits to the critical couplings; (left) plain 6th order polynomial, (center) ratio of a 4th to 2nd order 



polynomial and (right) physical fit according to Eq. (2.3). 



In addition we implementd a new fit strategy which consists in writing down the interpolating 
function in physical units ("physical" fit in Fig. 4 (right)) 

(2.2) 



r,(0)J l+AAt2/r/(At)+BMVr/(Ai)' 

with Tc = \/ {a{^)Lt) and A, B, C fit parameters. The implicit relation between j3c and /i^ is given 
by 

|2 /t2 1 D,,4 /t4 



a^(^c(M^))|2-loop =a^(^.(0))b-loopX ii^^^Jl^. 



(2.3) 



Also for the "physical" fit we get a ;^^/dof ^ 1. Moreover we have extrapolated the critical line 
down to r = axis, under the assumption that the physical fit gives the correct behavior of the 
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M/r^(0) 

Figure 5: The phase diagram obtained extrapolating the critical line down to T — axis by means of 



Eq. (2.2). 
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Figure 6: Extrapolation to real chemical potentials of the quadratic (with {ajj.)^ = (0.235/)^), sextic con- 
strained, ratio of polynomials and "physical" fits. 
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Figure 7: Comparison of our extrapolations with other determinations in the literature. For the sake of 
readability, our extrapolations have been plotted without error bands and labels, since they can be easily 
recovered from the previous figure. Legenda: D'Elia, Lombardo, Ref. [6,7]; Azcoiti et ai, Ref. [14]; Fodor, 
Katz, Ref. [15]; Kratochvila, de Forcrand, Ref. [16]. 
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critical line at real jj. down to T = (Fig. (^). We obtain the following estimate of the critical 
value of jj. on the T = axis: 

Ai = ^/c/B r,(0) = 2.5904(93) T,{0) . (2.4) 

We have so far discussed several interpolations of the data in the jJ.^ <0 region. A crucial 
question is to ascertain if the related extrapolations to /^^ > are consistent between them. Unfor- 
tunately, as one can inspect in Fig. 6, different interpolations lead to somewhat distinct extrapola- 
tions and, unless an extra-argument is found to make one fitting function preferable with respect 
to the others, one cannot rely on a unique extrapolation, except in the region ji/T < 0.6. In Fig. 7 
several determinations of the critical line existing in the literature are presented together our re- 
sults. Looking at Fig. 7, one could comment that the extrapolation of the "physical" fit exhibits the 
same trend as data from reweighting, whereas that from the sextic constrained fit mimics the strong 
coupling behavior [17], the other two extrapolations of ours lying in-between. However, previous 
determinations at real /i in the literature seem to be in fair agreement up to /^/r ~ 1 .2. If one takes 
this common trend as benchmark for our extrapolations, the "physical" and the polynomial ratio 
(4,2) seem to be favoured. 

We have tried to include in our fit also data at real chemical potential available from the lit- 
erature (see Fig. 7). A serious limitation of this combined approach is the inhomogeneity of the 
data presently available, due to different lattices and systematics. However, if the inhomogeneity of 
data at real pL will be reduced by new Monte Carlo determinations, the combined-fit strategy could 
bring along an appreciable improvement. 

3. Conclusions 

We have revisited the application of the method of analytic continuation from imaginary to 
real chemical potential in QCD with Nf = 4 degenerate flavors. The aims were: 

• to determine precisely the pseudo-critical hne Pdl^^) in the region of negative /i^ (20 data 
points almost uniformly distributed in the region — (7r/12)^ < (a/x)^ < 0); 

• to exploit interpolating functions sensitive to possible deviations of the critical line from the 
quadratic behavior in jj. for larger absolute values of /i (these deviations were clearly seen in 
QCD-like theories, such as 2-color QCD and finite isospin QCD, where it was given com- 
pelling evidence that their neglect could mislead the analytic continuation to real chemical 
potential); 

• to extrapolate the newly adopted interpolations to the region of real jj. and to re-determine, 
therefore, the critical line in QCD. 

We found that deviations from the quadratic behavior in jj. of j3c(/i^) at negative are visible 
in QCD with Nf = 4. Several kinds of functions able to interpolate them lead to extrapolations to 
real jU. which start diverging from each other for /i/T > 0.6. The shortcomings of the method of 
analytic continuation could be less severe for Nf = 2 or Nf = 2+1 (where sensitivity to nonlinear 
terms in /i^ could be enhanced). Moreover possible improvement could come by theoretical de- 
velopment able to discriminate between interpolations, or by a combined numerical strategy aimed 
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at gathering information from different approaches (such as reweighting, canonical approach, etc.) 
appUed so far independently from each other. 
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